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This paper deals with the problem of estimating the coupling constant S of a mixing quantum 
Markov chain. For a repeated measurement on the chain's output we show that the outcomes' time 
average has an asymptotically normal (Gaussian) distribution, and we give the explicit expressions 
of its mean and variance. In particular we obtain a simple estimator of 9 whose classical Fisher 
information can be optimized over different choices of measured observables. We then show that 
the quantum state of the output together with the system, is itself asymptotically Gaussian and 
compute its quantum Fisher information which sets an absolute bound to the estimation error. The 
classical and quantum Fisher informations are compared in a simple example. In the vicinity of 
S = we find that the quantum Fisher information has a quadratic rather than linear scaling in 
output size, and asymptotically the Fisher information is localised in the system, while the output 
is independent of the parameter. 

PACS numbers: 



I. INTRODUCTION 

Quantum Statistics started in the 70's with the dis- 
covery that notions of 'classical' statistics such as the 
Cramer-Rao inequahty, the Fisher information, have 
non-trivial quantum extensions which can be used to de- 
sign optimal measurements for quantum state estimation 
and discrimination [IHl]. Recently, statistical inference 
has become an indispensable tool in quantum engineering 
tasks such as state preparation [5l [6] , precision metrol- 
ogy [3 IH], quantum process tomography [51 [TU], state 
transfer and tclcportation jTTJ [T^] , continuous variables 
tomography [13, 14 . 

Quantum system identification (QSI) is a topic of par- 
ticular importance in quantum engineering and control 
[TSj where accurate knowledge of dynamical parame- 
ters is crucial. This paper addresses the QSI problem 
for Markov dynamics from the viewpoint of asymptotic 
statistics, complementing other recent investigations [TB]- 

[ig. 

We illustrate the concept of a quantum Markov chain 
through the example of an atom maser 10^: identically 
prepared d-level atoms (input) pass successively and at 
equal time intervals r through a cavity, interact with the 
cavity field, and exit in a perturbed state (output) which 
carries information about the interaction (see Figure pi) . 
Neglecting the internal dynamics of atoms and cavity, 
and taking the latter to be of dimension k < oo, the 
evolution can be described in discrete time and consists 
of applying an interaction unitary U G M{C^ (E) C*^) for 
each time interval r when a new atom passes through the 
cavity. If the incoming atoms are in the pure state ^ € C 
and the cavity is initially in some state ip g C*^, then at 



where f/^'-' is the copy of U acting on cavity and the atom 
which at time n is at position I on the right side of the 
cavity. 

Consider now that the interaction depends on some 
unknown parameter 9 such that U — Ue and correspond- 
ingly -0" = ipg. Our identification problem is to esti- 
mate 9, by measuring the output rather that the sys- 
tem (cavity) which may not be directly accessible. For 
simplicity we restrict ourselves to the case of estimat- 
ing one parameter (the interaction strength), such that 
Ue = exp{-i9H) where H € M(C'^ (g) C'') is a known 
hamiltonian, but the results can be extended to multiple 
parameters as well as continuous time [20] . The questions 



time n the output plus cavity state is ip" € 



^fj-v On 



V-" := C/(n) (V'®" ® ifi) := f/^^^ • • • C/^") (V-^" «) ifi) (1) 



Incoming atoms Optical cavity Outgoing atoms 

FIG. 1: Atom maser: identically prepared atoms pass suc- 
cessively through a cavity, interact with the cavity field, with 
the outgoing atoms carrying information on the dynamics. 



we want to address are: how much information about 9 
is contained in the output state ipg, and how can we 'ex- 
tract' it ? 

The standard approach to such questions goes via the 
quantum Cramer-Rao inequality which shows that the 
variance of any unbiased estimator obtained by measur- 
ing a copy of a state pg, is lower bounded by the inverse 
of the quantum Fisher information F{9)~^ ^^^. Al- 
though this bound is generally not attainable for a single 
copy, it is asymptotically attainable, i.e. there exist a 
sequence of measurements M„ on n identically prepared 



systems^ and estimators 0„ such that 

lim nE{en-ef = F{ey^. 

n— foo 

Unlike this case where the n-copies Fisher information 
scales linearly with the number of systems, the Fisher 
information of the correlated states "09 depends on the 
joint state rather than that of a single sub-system, and 
there is no straightforward argument to show that its 
rescaled version is asymptotically attainable. Moreover, 
in the case of multi-dimensional parameters, this ap- 
proach would run into the same problems as the inde- 
pendent copies model, for which it is well known that 
the Cramer-Rao bound is not attainable even asymptot- 
ically [T]. 

For these reasons, we will pursue an asymptotic anal- 
ysis based on the concept of local asymptotic normal- 
ity (LAN) |22] which was recently extended to quantum 
statistics ^23-26^ and used to solve the (asymptotically) 
optimal estimation problem for general multiparametric 
models. In this paper we extend quantum LAN from 
independent to finitely correlated quantum states, and 
in the same time generalise results on LAN [37] and the 
Central Limit Theorem (CLT) for classical Markov chains 
[28]. 

In Theorem|4]we show that, locally with respect to the 
parameter 0, the output state ^^ can be approximated 
by a one-parameter family of coherent states, cf. (15 1. 



This result can be seen as the Markov version of asymp- 
totic Gaussianity in coherent spin states (CSS) [25]. As 
a by-product, we obtain the asymptotic expression of the 
quantum Fisher information F = F{6) (per atom) of the 



output state (16), which is equal to the 'Markovian vari- 



ance' of the driving hamiltonian. This quantity should be 
understood as the limit of the rescaled quantum Fisher 
informations F^")(0) of the n atoms family of states f/'g 



no) 



lim 

n— 4-00 



i^(")(6') 



In particular, F{9)^^ sets an asymptotic lower bound on 
the variance nE(^„ — 6*)^ of any sequence of unbiased es- 
timators {9n}. Using the formalism of LAN we can show 
that the lower bound is attainable but at the moment we 
do not not know the explicit form of the optimal measure- 
ment, and we expect it to be non-separable, and possibly 
unfeasible with current technology. 

In Theorem |3] we analyse the more realistic set-up of 
repeated, separate measurements performed on the out- 
going atoms, and show that the mean of the outcomes 
is asymptotically Gaussian and can be used to estimate 
9. The corresponding classical Fisher information can 
be maximised over different measured observables, so 
that the experimenter can perform the most informative 
separable measurement, and compare its performance 
with the benchmark given by the quantum Fisher in- 
formation (see Example). This generalises Wiseman's 
adaptive phase estimation protocol where a particular 



field quadrature is most informative among all quadra- 
tures [30) . With the same techniques, similar results can 
be obtained for means of other functionals of the mea- 
surement data such as correlations between subsequent 
atoms. However it remains an open problem to find 
the (asymptotic) classical Fisher information contained 
in the complete measurement data. 

The next two sections introduce the key concepts un- 
derlying our results: local asymptotic normality and er- 
godicity. The main results are contained in Theorems 3 
and 4. To illustrate these results we analyse a simple ex- 
ample based on a XY interaction for which we compute 
the quantum Fisher information and we plot the classi- 
cal Fisher information for different output observables. 
An interesting feature of this model is the divergence of 
the asymptotic quantum Fisher information per atom at 
vanishing interaction, which is due to a quadratic rather 
than linear scaling of the 'usual' quantum Fisher informa- 
tion F'^"\9) for 6* « 0. We investigate this behaviour and 
find that with the scaling 9 — u/n, the model converges 
to a simple unitary rotation model on the system, with 
input passing into the output unperturbed. We conclude 
with a discussion on further extensions, open problems 
and connections with other topics. 



II. LOCAL ASYMPTOTIC NORMALITY IN 
CLASSICAL AND QUANTUM STATISTICS 

In this section we briefiy review some asymptotic 
statistics techniques, show how they extend to quantum 
statistics, and explain why this is useful. The aim is 
to introduce the concept of local asymptotic normality, 
which will be encountered in the main results. Theorems 
[3]and|l 



A. Asymptotic estimation in classical statistics 

A typical problem in statistics is the following: es- 
timate an unknown parameter 9 = {9i, . . . ,9p) G E^, 
given the random variables Xi , . . . , X„ which are in- 
dependent and identically distributed (i.i.d.) and have 
probability distribution Pg depending 'smoothly' on 9. 
If 9n ■— 6'„(Xi, . . . , X„) is a unbiased estimator, that is 
IE(0„) = 9, then the Cramer-Rao (C-R) inequality pro- 
vides the following lower bound to its (rescaled) covari- 
ance matrix 



nE 



9) >I{9)- 



(2) 



where I{9) is a p x p positive definite real matrix called 
the Fisher information matrix at 9 and quantifies the 
amount of 'statistical information' about 9 contained in a 
single sample from the distribution Pg. lipg{x) = d¥g/dfi 
denotes the density of Pg with respect to some reference 



measure /i, then the matrix elements of I{9) are given by 
Tta\ f ( .d\ogp0{x) d\ogpg{x) v fn t> \ 

li^Jtj — / Peix) 1^ Q^ ^J.[dx) ^Eg{eg^,ig^j) 

and depends only on the local behavior of the statistical 
model {Pg : 9 G W} around the point 9. To give a simple 
example, if Xi € {0, 1} are independent coin tosses with 
P6)[Xi = 1] = 0, then the mean 



1 " 

'^n — ^ / Xi 
r} i. — ^ 



is an unbiased estimator of 9 whose distribution is asymp- 
totically normal according to the Central Limit Theorem 

where A^(m, V^ is the normal distribution with mean m 
and variance V . A simple calculation shows that in this 
case I{9)~^ = 9{1 — 9) so that ^„ achieves the C-R lower 
bound. Interestingly, the Fisher information diverges at 
the boundary of the interval [0,1], but note that the 
bound is meaningful only for points in the interior of 
the parameter space. A similar situation will occur later 
in an example of a quantum Markov chain. 

While in general there might not exist any unbi- 
ased estimators achieving the C-R bound for a given n, 
the theory says that 'good' estimators 6'„(Xi, . . . , X„) 
(e.g. maximum likelihood under certain conditions) are 
asymptotically normal with 



fn{9n-9)^N{Q,I{9r^) 



(3) 



such that the C-R bound is asymptotically achieved [3T] . 

Le Cam went a step further and discovered a more fun- 
damental phenomenon called local asymptotic norm,ality 
(LAN), which roughly means that the underlying statis- 
tical model {¥g : 9 G W} can be linearised in the neigh- 
bourhood of any fixed parameter, and approximated by 
a simple Gaussian model with fixed covariance and un- 
known mean f32J. To explain this, let us first note that 
without loss of generality we can 'localise' 9, i.e. write it 
as 



where 6*0 can be chosen to be a rough estimator based 
on a small sub-sample of size n — ■n}~'^ with < e ^ 1, 
and u is an unknown 'local parameter'. By a simple con- 
centration of measure argument j25j one can show that 
with vanishing probability of error, the local parameter 
satisfies ||w|| < rf . For all practical purposes, we can 
then use the more convenient local parametrisation by 
u G MP and denote the original distribution P" , , ^ 
by Pn.u- Now, LAN is the statement that there exists 
randomisation (classical channels) r„ and Sn such that 



lim 
lim 

n— ^oo 



sup 

u\\ <n^ 

sup 

u\\ <n^ 



\Tn{Pn,u) -A/'(w)||i 
\Pn,u-S,,{Af{u))\\, 







where |1 ■ ||i denotes the Li-norm, and Pn,u and Af{u) 
are the probability densities of Vn,u and respectively 
N{u, I{9q)~^). Operationally this means that one can use 
the data Xi, . . . , X„ to simulate a normally distributed 
variable with density Af{u), and viceversa, with asymp- 
totically vanishing Li-error, without having access to the 
unknown parameter u. This type of convergence is strong 
enough to imply the previous results on asymptotic nor- 
mality and optimality of the maximum likelihood estima- 
tion, but can be used to make similar statements about 
other statistical decision problem concerning 9. We will 
now show that a similar phenomenon occurs in quan- 
tum statistics, and indicate how it can be used to finding 
asymptotically optimal state estimation protocols. 



B. The Cramer- Rao approach to state estimation 

Let us consider the problem of estimating a one- 
dimensional parameter 9, given n identical and indepen- 
dent copies of a quantum system prepared in the (pos- 
sibly mixed) state pg G M(C'^) depending smoothly on 
9. Following [THU [H] we analyse an unbiased estimator 
9n based on the outcome of an arbitrary measurement 
AI on the joint state pf". By the classical Cramer-Rao 
inequality, the mean square error (MSE) of 9n is lower 
bounded by the inverse (classical) Fisher information of 
the measurement outcome: 



E 



> hiiey 



Thus, a 'good' measurement is characterised by a large 
Fisher information, but how large can Im{0) be ? The 
answer is given by the notion of quantum Fisher infor- 
mation associated to the model {pg : 9 G M} which is 
defined as 

F{9):=TT{pgCl), 

where the symmetric logaritmic derivative (sld) Cg is the 
quantum analogue of ig and is the selfadjoint solution of 
the equation 



dpg 
d9 



{Cgpg + pgCg) 



As expected, for n identical copies the quantum Fisher 
information of the joint state {pf " : 6* e M} is nF{9) and 
the sld is given by 



£e(n) = ^£ 



(i) 



,(») 



with Cg acting on the i'th system. 

The Braunstein-Caves inequality shows that the quan- 
tum Fisher information is an upper bound to the classical 
one, i.e. Im{S) < nF{9) so that we obtain 



nE 



> F{oy 



Moreover, the constant on the right side can be 
achieved asymptotically by an adaptive measurement 
which amounts to measuring Cggin) for some point Oq 
which is a prehminary estimator of 9 obtained by mea- 
suring a small proportion of the systems. In summary, 
the optimal estimation rate for one dimensional param- 
eters is n~^F{6)~^, and can be achieved by means of 
separate measurements of the sld's Cg ■ 

Let us turn now to the case where pg depends on a 
multidimensional parameter 6 = (6*1, ... , 9p) g W. The 
quantum Fisher information matrix can be defined along 
similar lines and all previous inequalities hold as matrix 
inequalities, in particular the covariance matrix of an un- 
biased estimator 0„ satisfies 



lE 



-Of {On -6) >F{9) 



(4) 



However, unlike the classical case, and the one- 
dimensional quantum case, the right side is in gen- 
eral not achievable, even asymptotically ! This purely 
quantum phenomenon has a simple intuitive explana- 
tion: the optimal estimation of the different coordinates 
{6i, . . . ,9p) requires the simultaneous measurement of 
generally incompatible observables, the associated sld's 
(£0,i(n),. ..,£e,p(n)). 

Coming back to the original goal of estimating the pa- 
rameter 9, the above Fisher information analysis implies 
that the optimal measurement procedures must depend 
on the chosen figure of merit. Thus, one should not aim 
at saturating matrix inequalities such as Q but at find- 
ing asymptotically attainable lower bounds for the risk 
(multiplied by n) 



nE p„ - e)Gi9)i9n -9y^= nTr (^G(0)Var(0„) j , 

assuming for simplicity a quadratic loss function with 
positive weight matrix G{9). Taking the trace with G{9) 
on both sides of Q gives the generally non-attainable 
lower bound Tt{G(9)F{9)^^), and other examples can be 
derived from different versions of the quantum Cramer- 
Rao inequality such as Belavkin's right and left inequal- 
ities [5]. Holevo [T] derived a more general bound and 
showed that it is achievable for families of Gaussian states 
with unknown displacements, but until recently it re- 
mained an open question whether this bound was asymp- 
totically attainable for finite dimensional states. By fur- 
ther refining the techniques of the unbiased estimation 
set-up, Hayashi and Matsumoto "33] showed that the 
Holevo bound is indeed asymptotically attainable for gen- 
eral families of two-dimensional quantum states. Com- 
plementing this frequentist asymptotic analysis, Bagan 
and coworkers [34 solved the optimal qubit estimation 
problem for any given n in the Bayesian set-up with in- 
variant priors. However, neither of these approaches was 
successful in solving the (asymptotically) optimal estima- 
tion problem for general (mixed states) , multi-parametric 
models pe G M{C'') with fc > 2. 



C. Local asymptotic normality for quantum states 

At this point, a natural question to ask is whether the 
phenomenon of local asymptotic normality occurs also in 
quantum statistics, and whether it can be used to design 
asymptotically optimal measurement strategies. Recall 
that in the classical case, the main idea was that for 
large n, the i.i.d. model could be approximated by a 
Gaussian model, in the sense that each can be mapped 
approximately into the other by means of classical chan- 
nels. Building on earlier work by Hayashi [351 ES], the 
quantum version of LAN has been derived in a series of 
papers [S3H^ to which we refer for the details of the 
constructions. Here we only mention the general result, 
and discuss in more detail the special case of pure states 
models which is more relevant for the present work. 

As in the classical case, we can localise to a neigh- 
bourhood of 9o such that 9 ^ 9q + uj ^/n with ||m|| < n*^ 
for some small e > 0, and we denote pn.u '■= pf"l / ^- 
Then there exist channels (normalised, completely pos- 
itive linear maps) T„ and Sn between the appropriated 
spaces such that 

lim sup |lT„(p„,„) -Nu® *„|li = 0, 

lim sup ||p„,„-5'„(7V„(8)$„)||i =0 (6) 

"^°°||u||<n« 



(5) 



where {<&„ : u S M^} is a family of quantum Gaussian 
states of a continuos variables system, and \Nu : u G K^} 
is a family of classical Gaussian distributions. Moreover, 
each family has a fixed covariance matrix and the dis- 
placement is a linear transformation of the unknown pa- 
rameter u. With this tool at hand, one can prove that the 
Holevo bound is attainable through the following three 
steps procedure: first localise 0, then send the remain- 
ing states through the channel r„ and then apply the 
optimal measurement for the limit Gaussian model. 

We illustrate LAN for the simple case of a one dimen- 
sional family of pure states on C^ 



|V«):=exp(-z0J)|V'), 



(7) 



where J is a selfadjoint operator, which satisfies 
(i/'l Jj?/;) = 0. The (non-unique) sld is given by 

and the quantum Fisher information is proportional to 
the variance of the 'generator' J 



F(0) = (V'e|>C^|V'e)=4(V|J' 



F. 



(8) 



In the case of pure states, the limit Gaussian model 
consists of a family of coherent states so that LAN can 
be intuitively understood by analysing the intrinsic geo- 
metric structure of the quantum statistical model, which 
is encoded in the inner products between vectors corre- 
sponding to different parameters. Indeed, with the usual 



definition of the local states \tp. 
simple calculation shows that 



lim {^pn,u\i>n,v) ^ lim (V'|exp(i(u- w)/ynJ)|'0)" 
= exp(-(M - vfF/8) = {V2Fu\V2Fv) (9) 



where \u) is a coherent state of a one mode continuous 
variables (cv) system, with displacements {u\Q\u) — u 
and {u\P\u) — 0. This means that locally, the 'shape' of 
the statistical model for n states, converges to that of a 
family of coherent states. By using the central limit, we 
can identify the collective observables which converge "in 
distribution" to the coordinates of the cv system as 







V2nF 



C0„ (n) 



V2nF 



Ea 



(^) 



P 



Q 



such that the rescaled sld Csg{n)/y/n converges to the 
sld of the limit Gaussian model {\^/2Fu) : u e K}, as 
expected. With a more careful analysis of the speed of 
convergence in (lo]), it can be shown that the weak LAN 
can be upgraded to the strong version described by ([5| 
and ^ ^. 

The goal of this paper is to derive the weak LAN for 
the output state of a mixing Markov chain together with 
its classical counterpart for averages of simple measure- 
ments. The discussion around the i.i.d. models will hope- 
fully provide the necessary intuition about the statistical 
meaning of LAN in the Markov set-up and convince the 
reader that the quantity ( 16 1 plays the role of asymp- 
totic quantum Fisher information per atom. We leave 
the purely technical step for of deriving the strong LAN, 
and proving the achievability of the quantum Fisher in- 
formation for 1201. 



III. MIXING QUANTUM MARKOV CHAINS 

For later purposes we recall some ergodicity notions 
for a Markov chain with fixed unitary U. The re- 
duced n-steps dynamics of the cavity is given by the CP 
map p !—>■ T^{p) where T* is the 'transition operator' 
T*(p) = TvaiU \ip){ip\ <S) pW) with the trace taken over 
the atom. We say that T* is mixing if it has a unique 
stationary state T^,{pst) = pst, and any other state con- 
verges to Pst i.e. ||T'"(p) — Pst\\i — >■ 0. Mixing chains have 
a simple characterisation in terms of the eigenvalues of 
T^, , generalising the classical Perron- Frobenius Theorem. 

Theorem 1. T* is mixing if and only if it has a unique 
eigenvector with eigenvalue Ai = 1 and all other eigen- 
values satisfy \X\ < 1. 

As a corollary, the convergence to equilibrium is es- 
sentially exponentially fast \\T"-{p) — pst\\i < Cfcn'^|A2|" 



where A2 is the eigenvalue of T with the second largest 
absolute value [37]. The following theorem is a discrete 
time analogue of the perturbation theorem 5.13 of [55] 
and its proof is given in the appendix. 

Theorem 2. LetT{n) : Af(C*') ->: M{C'') be a sequence 
of linear contractions with asymptotic expansion 



1 



1, 



T{n) = To + -=Ti + -T2 + 0{n 
^n n 



-3/2N 



(10) 



such that Tq is a mixing CP map with stationary state 
Pst ■ Then Id — Tq is invertible on the orthogonal comple- 
ment of 1 with respect to the inner product {A, B)st '■= 
TT{pst-^ B). Assuming (l,Ti(l))st =Q we have 

lim T(n)"[l] =exp(A)-l, 

n— >oo 

where A = Tr (^.4(^2(1) + Ti o (Id - Tq)-^ o ri(l)) . 

From now on we assume that T = Tg is mixing. Since 
the cavity equilibrates exponentially fast we can choose 
the stationary state p^st ^'^ initial state, without affecting 
the asymptotic results below. 



IV. MAIN RESULTS 

This section contains the main results of the paper. 
The first subsection deals with estimators based on out- 
comes of separate identical measurements on the output 
atoms. In Theorem [3] we prove the asymptotic normal- 
ity of the outcomes' time average, and find the explicit 
expression of the asymptotic Fisher information of such 
statistics. Similar results can be obtained for time aver- 
ages of functionals depending on several outcomes, such 
as the correlations between subsequent atoms. In gen- 
eral these will provide higher Fisher information since 
the measurement outcomes are not independent but have 
exponentially decaying correlations. It remains an open 
problem is to find the 'full' Fisher information of the 
measurement stochastic process. The second subsection 
deals with the intrinsic statistical properties of the quan- 
tum model ipg . In Theorem [4] we prove that the quan- 
tum model is asymptotically normal, and we find the 
explicit expression of the quantum Fisher information F 
per atom. 



A. Simple measurements 

Consider a simple measurement scheme where an ob- 
servable A e M(C'') is measured on each of the outgoing 
atoms, and let A*^'^ be the random outcome of the mea- 
surement on the Vs atom. By stationarity, all expecta- 
tions ]E6/(A'^'') are equal to 

{A)e :- Tr (|V)(V'| ® pit ul{A (g> l)Ue 



and by ergodicity of the measurement process 



1 " 



(0 



{A)e, 



a.s. 



1=1 



Generically the right side depends smoothly on 9 and 
can be inverted (at least locally) so that 9 = /((A)e), 
for some well behaved function /, providing us with the 
estimator 9n '■= /(An)- As argued above, to analyse its 
asymptotic performance we can take 9 ~ 9q + u/\/n with 
9q fixed and u an unknown 'local parameter'. 

Theorem 3. Let 9 = 9o + uj %/n with 9q fixed and let 
A e M{<C'^)sa be such that (A)^,, = 0. Then ^JnAn is 
asymptotically normal, i.e. as n — > cx) 



E« 



i+tt/yn 



[exp((itV^A„)] -^ Fjt) 



^ ' — ^ is the character- 



where Fu{t) :— exp {i^{A)ut ^ 

istic function of the distribution N{^{A)u,a'^) with 

^l{A)■.^i{[H,A(g)l + l(^B])g„ (11) 

a^{A):={A^)g„+2{A®B)g, (12) 

B:={U-To)-'mulA^lUe,m (13) 

Note that for u — we obtain a quantum extension of 
the Central Limit Theorem for Markov chains j28| . 

Before proving the theorem we show that 0„ is asymp- 
totically normal and find its mean square error. By ex- 
panding 9n '■= ./(An) around {A)g„ and using the prop- 
erty f'{{A)g) = {d{A)g/d9)-^ we have 

0n = fi{A)e,) + /'((A)eJ(A„ - {A)e„) + Or{n-') 
^9o + fi{A)-^A„ + Op{n-^) 



so that 



lim nK 



= hm E 



y/nA-n 



a\A) 



The limit can be seen as the inverse (classical) Fisher 
information per measured atom, which in principle can 
be minimized by varying A and/or the input state '0. 



Proof. By using ([T]) we can rewrite the characteristic 
function as 



Ee [exp((zt^A„)] = (exp(iiV^A„))g = (r(7i)"[l])) 
where T{n) : M{€^) -^ M{€^) is the map 
T{n) : X ^-> Eso \e"''/^' L^tA/^ ^ ^\ ^-iHu/^i 



with Effg [Y\s] := (i/'If^J ^^SolV') the 'conditional expec- 
tation' onto the system. We expand T{n) as in ( 10 ) with 



To{X):^Ee,[l<E>X\s] 

Ti{X) ■.^iEgJu[H,l(g)X]+tA(g,X\s] 

u^ 

T^{X):^--Eg,{{H\U,\®X\M 

+2 

--Ee„ {A^®X\^\ -utEg^\\H,A®X\\^\ 

Since Tr(pstTi(l)) = {A)g^ = we can apply Theorem[2] 
so we only need to compute the coefficient A. The first 
part is Tr(p,tr2(l)) - ~e{^A^)gj2 - ut{[H,A® l])g, 
and the second part is iiTr(psi • Ti{B)), with B defined 
in ([13]). 

n 



B. Quantum Fisher information and LAN 

Recall that the joint state of the system and output 
atoms is 

^^ = Ug{n) (V-®" ® v) := C/(i) • • • [/(") (V-®" ® ^) 
= exp(-ie'ij(i))---exp(-i6li/(")) {ip®'' ® if) 

where -/?*■*' represents the copy of the interaction hamil- 
tonian which acts on the system and the i'th atom, and 
is ampliated by the identity on the rest of the atoms. 
It is important to note that in general the commutants 
[H^''\ H^J")] are nonzero for i ^ j since both hamiltonians 
contain system operators. This means that the model -tjjg 
is not a covariant one as ([7]), i.e. we cannot write it as 
-0g = exp{—i9H{n)) {tp^^ ® if) for some 'total hamilto- 
nian' H{n). In particular, as we will see, the quantum 
Fisher information depends on 9. However, we can write 



dpl 
d9 



niAY ■ where p^ = IV-g ) (V'e I and 



Hg{n)^Y.Hr ^-T.ulr ■■■Ut''H^''U, 



[He{n),p'^] 



(1) ...r/('"^'77"(*)r/''"^'* •••r/^^) 



As in ([8]) it follows that the quantum Fisher information 
F(") (0)is 

F^-H9)^A{Hg{nf)g^^-A{Hg{n))l^. 

Here and in the next few lines the expectation {■)g.n is 
taken with respect to the state tp'g, rather than the sta- 
tionary state whose expectation is denoted {■)g. We will 
see that in asymptotics, we can revert to the stationary 
state expectation. Ignoring for the moment the second 



term on the right side, we write 



(He{nr)e^n ^ E ( (^< 



0,n 



+ E ({^i^^i^-^}) 



E 

l<i<j<n 



unitary famihes of pure states, as variance of the driv- 
ing hamiltonian as shown in section III] By extending the 
weak convergence to strong convergence as described in 
section ITT] and applying similar techniques as in the i.i.d. 
case pS] it can be shown that there exists a two step 
adaptive measurement which asymptotically achieves the 
smallest possible variance 



Now, by ergodicity we have 



lE 



P- 



9„ := 9o + Un/Vn. (17) 



lini — 



Z^ 



H. 



ii) 



{H')e 



Similarly, by using the Markov property, it can be shown 
that 



lim - y 

l<'i<^'<n 






i=0 



Since for coherent models the optimal measurement is 
that of the canonical variable which is conjugate to the 
'driving' one, we conjecture that quantum Fisher infor- 
mation is achieved by a measuring a collective observable 
which is conjugate to the hamiltonian iJ in a more gen- 
eral Markov CLT than that of Theorem [S] 

Proof of Theorem 4- By (JT]) we can rewrite the inner 
product as 



where K := {ip\H\ij) € M{C''). Assuming that {H)g = 
0, the bias term {He{n))g ^ is sub-linear and from the 
above limits we obtain 



iVo 



/sl^, 



eo+«/Vnl^So+f/V"' 



= Tr(pr(n)"[l]) 



where 



1 



lim -P^''>{eo) = Pi9o) 



(14) 



T{n) -.Xy-^Ee 



^ff/v^(l ^ x)e"™-^/^ 



with P{9o) the asymptotic quantum Fisher information Us expansion is 



per atom, given by (16). The next theorem strengthens 
this conclusion, by showing that the quantum model ipg 
converges (locally around any Oq 7^ 0) to a coherent state 
model with quantum Fisher information P{0q). 

Theorem 4. Suppose that {H)g^ = 0. Then the family 
of (local) output states {ip' 



To{X):=Eg,[l®X\s\ 

Ti{X) -.^iEg, [uH{l ®X)~ v{l ® X)H\s] 

T2{X) := ^' 



Oq+u/^ 



u £ M} converges 
to a family of coherent states {\^P/2u) : u G M} with 
Pd displacement {^/F/2u\Q\^/fJ2u) = ^/F/2u. More 
precisely, as n ^>- 00 



-Eg, [u^H^il (g,X)+ v^{l (g, X)H^\s\ 
-Eg, [uvH{l(E)X)H\s]. 



W 



o+u/^\'^eo+v/^) 



^ia{u —V ) J 



^ gia(«^-t,^)g-F(u-«)V8 



(15) 



where a £M. is a constant and P = P{Oq) is the asym,p- 
totic quantum Pisher information (per atom) 

Pi9o) := 4 [{H^)g„ + {{H, (Id - To)-\K)})g,] , (16) 

with K := Eg„[H\s] = {il:\H\%l)) . 

Note that e*"'^" ^" > is an irrelevant phase factor which 
can be absorbed in the definition of the limit states. The 
claim that P is the asymptotic quantum Fisher informa- 
tion (per atom) of the output follows from the conver- 
gence (151 to the one dimensional coherent state model 
and the fact that the latter has quan- 



{WPl2u) :wg 

turn Fisher information F, cf. [HIH] , and agrees with the 

limit (14). Note also that similarly to (12), the left side 



of ( 16 ) can be identified with the asymptotic variance 



appearing in the CLT for the operator 2H . This agrees 
with the formula of the quantum Fisher information for 



and the condition of Theorem 2 holds: Tr(psfTi(l)) — 
{H)g„ = 0. Since Ti(l) = i{u — v)K, the contribution 
from Ti^pstTi o (Id — Tq)~^ o Ti(l)) can be written as 

-{u- vf Re(iJ(Id - To)-'{K))g, 
- t{u^ - v^) Im(iJ(Id - To)-\K))g„ 

Moreover, Tr(p,tT2(l)) = -{u- v)^H^)gj2, so adding 
the two contributions we obtain the desired result. 

n 



V. EXAMPLE 

In this section we illustrate the main results on the 
classical and quantum Fisher informations, for a simple 
discrete time XF-interaction model. We also analyse 
the behaviour of the quantum Fisher information in the 
vicinity of the point 6* = at which the chain is not 
ergodic, and find that it is quadratic rather than linear 
in n and is concentrated in the system rather than the 
output. 

Let Ug — exp(—i6H) with H the 2-spin 'creation- 
annihilation' hamiltonian H = i(cr+ (g) cr_ — cr_ ® (7+) 




o.oW' 

0.0 



FIG. 2: Classical Fisher information for spin measurements 
as function of spin components Uy , riz 



where a± are the raising and lowering operators on C^, 
and let \tp) := a\0) + 5e*^|l) be the input state, with 
a^ + b^ = 1. The transition matrix Tq has eigenvalues 
{1, c, {c{c + 1) ± v/c2(l - c)2 - 16a262c(l - c2))/2}, and 
is mixing if and only if c := cosSq 7^ 1- The quantum 
Fisher information (|16|) is 



F := 



Wa'^b'^ 



(l-c)(l-c + 4a262c) 



(18) 



which is independent of the phase /. This can be com- 
pared with the value of the classical Fisher information 
^{(jfi)'^ / a'^ {af{) obtained by measuring the spin in the di- 
rection n for each outgoing atom, cf. Theorem ^ For 
c — 0.5, / = 0, 6 = 0.8 the quantum Fisher information is 
F = 5.03 while that of the spin measurement varies from 
to 1.15 (see Figurel2]). 



An interesting, and perhaps surprising feature of the 
quantum Fisher information ( 18 ) is that it diverges 



at vanishing coupling constant, due to the factor 1 — 
cos{9) ss 02^2 in the denominator. This singularity 
arises from the second term on the right side of ( 16 1 and 
stems from the fact that the Markov chain is not mix- 
ing at 9 — 0. To get an intuition for this phenomenon 
let us compute (according to the standard methodology 
[U El [21] ) the quantum Fisher information at 9 = for 
the family 

where the input state is chosen to be |'0) = (|0) + |l))/\/2 
and the initial state of the 'cavity' is \(p) = |0). Since 



-iH{n) 



dUe{n) 



d9 



^ ( cr+ «) (t'l' - cr_ (g) cr 



.(') 



9=0 i=i 

one can easily verify that 

(Ve%|ff(n)|Ve=o>=0, 



and hence the quantum Fisher information at = is 

F^"^ = 4(V^,%|i7(n)2|V;,%) ^ 4||i/(n)|^,%)||2. 

Since C-lO) = the latter reduces to computing the 
squared norm of the vector 



y^ c{ii,- ■.i„)\ii,.-.in)- 



where ij € {0, 1} are the basis indices. Since a basis 
vector containing p indices equal to 0, can be obtained 
in p different ways by applying the lowering operator to 
an input tensor, the coefficients are 

n 

c(ii,...,*„) = (n-^z,)/2"/2. 




and the Fisher information is 



pi") 



-"^p2 
p=0 



n{n + 1) 



Thus, in contrast to the case of independent systems, the 
Fisher information scales quadratically rather than lin- 
early with the number of systems, hence the divergence 
of the quantum Fisher information per atom which repre- 
sents the asymptotic value of F^"'> /n. Note that similar 
quadratic scaling of the quantum Fisher information is 
encountered in phase estimation [7j and more generally 
in optimal estimation of unitary channels |40j . with the 
difference that it holds for any parameter, rather than at 
a single point. 

We take a closer look at the the states tpg by scaling 
the parameter as 9 = u/n as suggested by the Fisher 
information. This means that we know 9 with an accu- 
racy of n~^ (rather than the usual n~^''^) and we would 
like to find if the the estimation of u 'stabilises' in the 
asymptotic regime. By using the same technique as in 
Theorem |4] we have 



(C/nlC/J = (0|T(n)"[l]|0) 
where T{n) : M{C^) -^ M(C^) is the map 

T{n) : X ^ (V'|e™^/"(1 X)e-™^/"| 
We expand T{n) as 



(19) 



T{n) = Id 



Tn 



+ o{n 



where 



Tq: X^i{uKX -vXK), K := {^\H\i>). 



Plugging into (19) we obtain the limit 

lini(C/„IC/J = (0|exp(To)[l]|0) 

= {0\exp{i{u-v)K)\0). (20) 
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This agrees with the quadratic scahng of the quantum 
Fisher information at 6 = and shows that the quan- 
tum statistical model has a limit provided that the right 
scaling of parameters is used. 

Let us now look at the n-steps reduced dynamics of the 
system for the same scaling 6 = u/n with initial state |0). 
After a similar computation (with u — v) we obtain 

lim p{n) = lim T,(n)"(|0)(0|) 

n— J-oo n— >-oo 

= cxp{-iuK)\0){0\exp{iuK), (21) 

which means that for large n the system is effectively 
unitarily driven by the 'hamiltonian' K^ even though we 
started with an open system dynamics! This effect is in- 
teresting in itself and is reminiscent of the quantum Zeno 
effect. From (20) and (21) we can conclude that asymp- 



totically, the system and the output have pure states, and 
moreover the input passes undisturbed into the output 

IC/„)«IV')^"^exp(z^X)|0). 

In conclusion, unlike the ergodic case where the output 
contains information about the parameter, it is the sys- 
tem's state which carries all the information, and suc- 
cessive time steps amount to a simple unitary rotation. 
In particular, this explains the quadratic scaling of the 
quantum Fisher information F^"^' . In conclusion, the 
non-ergodic set-up exhibits interesting statistical features 
and should be analysed on its own in more detail. In 
particular, for practical applications it is important to 
see whether the quadratically scaled Fisher information 
is achievable asymptotically. 



VI. CONCLUSIONS AND OUTLOOK 



When the chain is not-ergodic, the Fisher information 
may exhibit a qualitatively different behaviour, such as 
quadratic rather than linear scaling with the number of 
output systems, which bears some similarity with that 
encountered in phase estimation [7 . Understanding this 
behaviour in a more general scenario, and the possible 
applications in precision metrology are topics for future 
research. 
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Appendix: Proof of Theorem [2] 

Since Tq is a mixing CP-map, the identity is the unique 
eigenvector with eigenvalue 1 and all other eigenvalues 
have absolute values strictly smaller than 1. For n 
large enough T{n) has the same spectral gap property 
and we denote by (A(n),a:(n)) its largest eigenvalue and 
the corresponding eigenvector such that \{n) — > 1 and 
x{n) -^ 1. 

Since T{n) is a contraction 



|r(n)"(l)-r(n)"(x(n))|| < ||r(n)"|| ||l-a;(n)|| 



0. 



On the other hand since T{n)'^{x{n)) = X{n)^x{n), we 
have T(n)"(l) ~> exp(A)l, provided that the following 
limit exists 

lim A(n)" := exp(A). 



Quantum system identification is an area of signifi- 
cant practical relevance with interesting statistical prob- 
lems going beyond the state estimation framework. We 
showed that in the Markovian set-up this problem is very 
tractable thanks to the asymptotic normality satisfied 
by the output state and the time average of the sim- 
ple measurement process. This may come as a surprise 
considering that the output is correlated, but is in per- 
fect agreement with the classical theory of Markov chains 
where similar results hold [27 . The theorems can be ex- 
tended to strong LAN with multiple parameters, contin- 
uos time dynamics, and measurements on several atoms 
[20j . However, as in process tomography, full identifica- 
tion of the unitary requires the preparation of different 
input states. The optimisation of these states, and the 
case of non-mixing chains are interesting open problems. 
Another open problem is to find the classical Fisher infor- 
mation of the simple measurement process, rather than 
that of time averaged functionals. Since Markov chains 
are closely related to matrix product states [HI |42] , our 
results are also relevant for estimating matrix product 
states ESI. 



We prove that this is the case by using the Taylor ex- 
pasions 

A(n) = Ao + ^Ai + -A2 + 0(n-3/'), 
y/n n 



1 



x{n) = xo + ^xi + -X2 + 0{n •^'^). 
y/n n 

Then we can solve the eigenvalue problem in successive 
orders of approximation 

To{xo) = Ao^o, 
Tq{xi) +Ti{xo) = AqXi+AiXo, 
Toix2) +Ti{xi) +T2{xo) = Aoa;2 + A1.T1 + A2a;o (*). 

From the first equation we have Aq = 1, 2:0 = 1- Inserting 
into the second equation we get 



(To-Id)(xi 



-(Ti-Aild)(l), 



and by taking inner product with 1 we obtain Ai = 
(l,Ti(l))st — 0, by using the assumption. Hence 



Xi 



iId~Tor'oTi{l) + cl 



10 



where (Id — Tq)~^ denotes the inverse of the restriction 
of Id — To to the orthogonal complement of 1, and c is a 
constant. 

Similarly, from the third equation in (*) 

(To - Id)(x2) = -Ti{xi) + (Aald - T2)(l). 

which implies 

A2 = (l,T2(l)),t + (l,Ti(xi)),t 

= (1, T2(l) + Ti o (Id - To)-i o ri(l)),t 



Finally, 



lim A(n)" = lim ( Aq + 4^ + — + 0(n-3/2) 
= lim ( 1 H = cxp(A2). 

ri->-oo \ n J 



n 



[1] A. S. Holevo, Probabilistic and Statistical Aspects of 

Quantum Theory (North-Holland, 1982). [23 

[2] C. W. Helstrom, Quantum Detection and Estimation [24 

Theory (Academic Press, New York, 1976). 
[3] V. P. Belavkin, Theor. Math. Phys. 26, 213 (1976). [25 

[4] H. P. Yuen and Lax, M., IEEE Trans. Inform. Theory 

19, 740 (1973). [26 

[5] H. Haffner, W. Hansel, C. F. Roos, J. Benhelm, D. Chek- 

al kar, M. Chwalla, T. Korber, U. D. Rapol, M. Riebe, [27 

P. O. Schmidt, et al.. Nature 438, 643 (2005). 
[6] G. Breitenbach, S. Schiller, and J. Mlynek, Nature 387, [28 

471 (1997). 
[7] V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, [29 

1330 (2004). [30 

[8] B. L. Higgins, D. W. Berry, S. D. Bartlett, H. M. Wise- [31 

man, and G. J. Pryde, Nature 450, 393 (2007). 
[9] M. Lobino, C. Kupchak, E. Figueroa, and A. I. Lvovsky, [32 

Phys. Rev. Lett. 102, 203601 (2009). 
[10] M. Brune, J. Bernu, C. Guerlin, S. Deleglise, C. Sayrin, [33 

S. Gleyzes, S. Kuhr, I. Dotsenko, J. M. Raimond, and [34 

S. Haroche, Phys. Rev. Lett. 101, 240402 (2008). 
[11] B. Julsgaard, J. Sherson, J. I. Cirac, J. Fiurasek, and [35 

E. S. Polzik, Nature 432, 482 (2004). 

[12] J. F. Sherson, H. Krauter, R. K. Olsson, B. Julsgaard, 

K. Hammerer, J. I. Cirac, and E. S. Polzik, Nature 443, [36 

557 (2006). [37 

[13] K. Vogel and Risken, H., Phys. Rev. A 40, 2847 (1989). 

[14] L. G. Lutterbach and L. Davidovich, Phys. Rev. Lett. [38 
78, 2547 (1997). 

[15] H. Mabuchi and N. Khaneja, International Journal of [39 
Robust and Nonhnear Control 15, 647 (2005). 

[16] H. Mabuchi, Quatum Semiclass. Opt. 8, 1103 (1996). 

[17] M. Howard, J. Twamley, C. Wittmann, T. Gaebel, [40 

F. Jelezko, and J. Wrachtrup, New J. Phys. 8, 33 (2006). [41 
[18] S. G. Schirmer and D. K. L. Oi, Laser Physics 20, 1203 

(2010). [42 

[19] D. Burgarth and K. Maruyama, New J. Phys. 11, 103019 

(2009). [43 

[20] M. Guta and L. Bouten, in preparation. 
[21] S. L. Braunstein and Caves C. M., Phys. Rev. Lett. 72, 

3439 (1994). 
[22] A. van der Vaart, Asymptotic Statistics (Cambridge Uni- 



versity Press, 1998). 

M. Guta and J. Kahn, Phys. Rev. A 73, 052108 (2006). 
M. Guta and A. Jengova, Commun. Math. Phys. 276, 
341 (2007). 

M. Guta, B. Janssens, and J. Kahn, Commun. Math. 
Phys. 2V7, 127 (2008). 

J. Kahn and M. Guta, Commun. Math. Phys. 289, 597 
(2009). 

R. Hopfner, J. Jacod, and L. Ladelli, Probab. Theory 
Related Fields 86, 105 (1990). 

S. Meyn and R. L. Tweedie, Markov chams and Stochas- 
tic Stability (Cambridge University Press, 2009). 
J. M. Radcliffe, J. Phys. A: Gen. Phys. 4, 313 (1971). 
H. M. Wiseman, Phys. Rev. Lett. 75, 4587 (1995). 
E. L. Lehmann and G. Casella, Theory of point estima- 
tion (Springer, 1998). 

L. Le Cam, Asymptotic Methods in Statistical Decision 
Theory (Springer Verlag, New York, 1986). 
M. Hayashi and Matsumoto, K., quant-ph/0411073. 
E. Bagan, Ballester, M. A., GiU, R. D., Monras, A., and 
Munoz-Tapia, R., Phys. Rev. A 73, 032301 (2006). 
M. Hayashi, Bulletin of the Mathematical Society of 
Japan 55, 368 (2003), ( in Japanese; Translated into En- 
glish in | quant-ph/0608198 ) . 
M. Hayashi. 

B. M. Terhal and D. D. P., Phys. Rev. A 61, 022301 
(2000). 

E. Davies, One-parameter semigroups (Academic Press 
Inc (London) Ltd, 1980). 

B. Kiimmerer and H. Maassen, Infinite Dimensional 
Analysis, Quantum Probability and Related Topics 3, 
161 (2000). 

J. Kahn, Phys. Rev. A 75, 022326 (2007). 
M. Fannes, B. Nachtergaele, , and R. F. Werner, Com- 
mun. Math. Phys. 144, 443 (1992). 

C. Schon, E. Solano, F. Verstraete, J. I. Cirac, and 
M. Wolf, Phys. Rev. Lett. 95, 110503 (2005). 

M. Cramer, M. B. Plenio, S. T. Flammia, R. Somma, 

D. Gross, S. D. Bartlett, O. Landon-Cardinal, D. Poulin, 
and Y.-K. Liu, Nature Commun. 1, 149 (2010). 



